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ABSTRACT 

The simultaneous detection of electromagnetic and gravitational wave emission from merging neu- 
tron star binaries would aid greatly in their discovery and interpretation. By studying turbulent 
amplification of magnetic fields in local high-resolution simulations of neutron star merger conditions, 
we demonstrate that magnetar- level (> 10 16 G) fields are present throughout the merger duration. 
We find that the small-scale turbulent dynamo converts 60% of the randomized kinetic energy into 
magnetic fields on a merger time scale. Since turbulent magnetic energy dissipates through reconnec- 
tion events which accelerate relativistic electrons, turbulence may facilitate the conversion of orbital 
kinetic energy into radiation. If 1CU 4 of the ~ 10 53 erg of orbital kinetic available gets processed 
through reconnection, and creates radiation in the 15-150 keV band, then the fluence at 200 Mpc 
would be 10 _7 erg/cm 2 , potentially rendering most merging neutron stars in the advanced LIGO and 
Virgo detection volumes detectable by Swift BAT. 

Subject headings: magnetohydrodynamics — turbulence — stars: neutron — magnetic fields - 
hydrodynamics — gamma-rays: bursts — X-rays: general — gravitational waves 



1. INTRODUCTION 

The in-spiral and coalescence of binary neutron star 
systems is a topic of increasingly intensive research in 
observational and theoretical astrophysics. It is antic- 
ipated that the first direct detections of gravitational 
wave (GW) will be from compact binary mergers. Binary 
neutron star (BNS) mergers are also thou ght to produce 
short - hard gamma-ray bursts (SGRB 's) (Eic hler et al. 
19891 INarayan et al|[l99"2"l |Belczynski et al.||2006| |Met- 
zger et al.||2008p . Simultaneous detections of a prompt 
gravitational wave signal with a spatially coincident elec- 
tromagnetic (EM) counterpart dramatically increases the 
potential science return of the discovery. For this rea- 
son, there has been considerable interest as to which, if 
any, detectable EM signature may result fro m the merger 
2012] |Piran et al.||2013[ ). Other than 



(Metzger & Berger 



SGRBs and their a: 



tcrglows, including those viewed off- 



axis ( Rhoads|1997 van Eerten fc MacFadyen|20TT ), sug- 
gestions include optical afterglows associated with the 
rad io-active decay of tidally expelled r-proce ss mate- 
rial( |Li fc Paczynski|[l998| |Metzger etld1j20iq | ) (though 
detailed calc ulations indicate they are taint ( Barnes fc 
Kasen 2013)), radio afterglows following the interaction 
of a mildly relativistic shell with the interstellar medium 
(Nakar & Piran [2TjlT| ), and high -energy pre-merger emis- 
sion from resistive magnetos phere interactions (Hansen 
fc Lyutikov|[200T| |Piroj|2012| ) 



Merging neutron stars possess abundant orbital ki- 
netic energy (~ 10 53 ergs). A fraction of this en- 
ergy is certain to be channelled through a turbulent 
cascade triggered by hydrodynamical instabilities dur- 
ing merger. Turbulence is known to amplify magnetic 
fields by stretching and folding embedded field lines in 
a process known as the small-scale turbulent dynamo 



( Vamshtci nfe Zel'dovich|1972||Ch ertkov et al.|1999| [To 



bias et al. 2011 " Beresnyak 2012 H Amplification stops 
when the magnetic energy grows to equi partition with 
the energy containing turbulent eddies ( |Schekochihin 



et al.|2002[|2004||Federrath et al.|20TT| ). An order of mag- 
nitude estimate of the magnetic energy available at sat- 
uration of the dynamo can be informed by global merger 
simulations. These studies indicate the presence of tur- 
bulence following the nonlinear saturation of the Kelvin- 
Helmholtz (KH) i nstability activated by shearing at the 
NS surface layers ([Price fc Rosswog|2006[ |Liu et al.|2008 " 



Anderson et al.|[20U8| |Rezzolla et al. 2011 " Hiacomazzo 
et al.||201ip . The largest eddies produced are on the ~ 1 
km scale and rotate at ~ 0.1c, setting the cascade time 
Wdy and kinetic energy injection rate Sk = -Eft" /Wdy 
at 0.1ms and 5 x 10 50 erg/s respectively. When kinetic 
equipartition is reached, each turbulent eddy contains 
5 x 10 46 erg of magnetic energy, and a mean magnetic 
field strength 



B 



[(MS 



> 10 16 G 



10 13 g/cm 3 



1/2 



/^eddy\ 

V 0.1c / 



(1) 



Whether such conditions are realized in merging neu- 
tron star systems depends upon the dynamo saturation 
time i sa t and equipartition level E^jEn- In particu- 
lar, if t sat < i mer ge then turbulent volumes of neutron 
star material will contain magnetar-level fields through- 
out the early merger phase. Once saturation is reached, a 
substantial fraction of the injected kinetic ener gy, 0.7ek , 



is resistively dissipated (Haugen et al. 2004) at small 



scales. Magnetic energy dissipated by reconnection in op 
tically thin surface layers will accelerate relativistic elec 
trons ( Blandford fc Eichler||1987[ Lyutikov & Uzdensky 
2003 1 , potentially yielding an observable electromagnetic 



counterpart, independently of whether the merger even- 
tually forms a relativistic outflow capable of powering a 
short gamma-ray burst. 

In this Letter we demonstrate that the small-scale tur- 
bulent dynamo saturates quickly, on a time t sat ^ Emerge, 
and that > 10 16 G magnetic fields are present through- 
out the early merger phase. This implies that the mag- 
netic energy budget of merging binary neutron stars is 
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Figure 1. Time development of volume-averaged kinetic and 
magnetic energies at resolutions between 16 3 and 1024 3 . Lower 
resolutions are shown in red and graduate to black with higher 
resolution. Top: The root mean square magnetic field strength in 
units of 10 16 G. When a turbulent volume is resolved by 16 3 zones, 
the small-scale dynamo proceeds so slowly that almost no amplifi- 
cation is observed in the first 1ms. Middle: The magnetic energy 
in units of the rest mass p c 2 shown on logarithmic axes. It is clear 
that the linear growth rate increases at each resolution. Bottom: 
The kinetic energy (upper curves) shown again the magnetic en- 
ergy (lower curves) again in units of p c 2 . For all resolutions, the 
kinetic energy saturates in less than 1 i c ddy 



controlled by the rate with which hydrodynamical insta- 
bilities randomize the orbital kinetic energy. Our results 
are derived from simulations of the small scale turbulent 
dynamo operating in the high-density, trans-relativistic, 
and highly conductive material present in merging neu- 
tron stars. We have carefully examined the approach to 
numerical convergence and report grid resolution crite- 
ria sufficient to resolve aspects of the small-scale dynamo. 
Our Letter is organized as follows. The numerical setup 
is briefly described in Section 2. Section 3 reports the 
resolution criterion for numerical convergence of the dy- 
namo completion time and the saturated field strength. 
In Section 4 we asses the possibility that magnetic recon- 
nection events may convert a sufficiently large fraction of 
the magnetic energy into high energy photons to yield a 
prompt electromagnetic counterpart detectable by high 
energy observatories including Swift and Fermi. 



2. NUMERICAL SCHEME AND PHYSICAL SETUP 

The equations of ideal relativistic magnetohydrody- 
namics (RMHD) have been solved on the periodic unit 
cube with resolutions between 16 3 and 1024 3 . 



V U N U = 
dB 



dt 



= V x (v x B) 



(2a) 
(2b) 

(2c) 



Here, 6^ = F^u" is the magnetic field four-vector, and 
h* = 1 + e* + p* / p is the total specific enthalpy, where 
b 2 /2 is the total pressure, p g is the gas pres- 



P 



Pg 



sure and e* = e + b~/2p is the specific internal en- 
ergy. The source term = pa^ — p(T/T ) 4 u^ includes 
injection of energy and momentum at the large scales 
and the subtraction of internal energy (with parameter 
T = 40MeV) to permit stationary evolution. Vortical 
modes at k/2ir < 3 are forced by the four- acceleration 

field a M = ^j- which smoothly deco rrelates over a large- 
eddy turnov er time, as described in |Zrake fc MacFadyen| 
( [2012] [2013] ). 

We have employed a realistic micro-physical equation 
of state (EOS) appropriate for the conditions of merg- 
ing neutron stars. It includes contributions from high- 
density nucleons according a relativist ic mean field model 



(Shen et al. 1998 Shen et al. 2010), a relativistic and 



degenerate electron-positron component, neutrino and 
anti-neutrino pairs in beta equilibrium with the nucle- 
ons, and radiation pressure. For our conditions, all the 
components make comparable contributions to the pres- 
sure. We have also employed a simpler gamma-law EOS 
and found close agreement for the conditions explored in 
this paper, indicating that the MHD effects are insen- 
sitive the EOS for trans-sonic conditions. The models 
presented in our resolution study use the far less expen- 
sive gamma-law equation of state. 

All of the simulations presented in this study use 
the H LLD approximate Riemann solver (Mignone et al. 
2009[ ) , which has been demonstrated as crucial in provid- 
ing the correct spa tial distribution of magne tic energy in 
MHD turbulence (iBeckwith & Stone||2011|). The solu- 



tion is advanced with an unsplit, second-order MUSCL- 
Hancock scheme. Spatial reconstruction is accomplished 
with the piecewise linear method configured to yield the 
smallest possible degree of numerical dissipation. The di- 
vergence constraint on the magnetic field is maintained to 
machine preci sion at cell c orners using the finite volume 
(| 2002 ). Full details of the numerical 
TnTZrake & MacFadyenl (|2012h. 



CT method of |T6th 
scheme may be foun 



3. RESULTS 



3.1. Growth and saturation of the magnetic energy 

Magnetic fields are amplified in our simulations by the 
small-scale turbulent dynamo. Turbulent fluid motions 
stretch and fold the magnetic field lines , causing expo - 



nential growth of magnetic energy (e.g. Moffatt 



1978). 



This growth is attributed to the advection and ditfu- 
sion of B through the MHD induction equation (Eq. 
[2b). When the magnetic field is weak (Em -C Ek) B 
evolves passively, and the turbulence is hydrodynamical. 
This limit is referred to as small-scale kinematic turbu- 
lent dynamo, and is well described by Kazentzev's model 
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Figure 2. Time development of the spectrum of kinetic (left) and magnetic (right) energy, given in erg/cm 3 /m — 1 for resolution 512 3 . 
During the startup phase the hydrodynamic cascade is not yet fully developed, this phase lasts for lc ddy The kinematic phase refers 
to the time when the hydrodynamic cascade is fully developed, but the magnetic energy is still energetically sub-dominant. During the 
kinematic phase, the the kinetic energy power spectrum is consistent with a Kolmogorov cascade, having a slope k~ 5 ' 3 over intermediate 
wavelengths. Meanwhile, the magnetic energy spectrum is consistent with the Kazentsev description, having a positive slope fc 3 / 2 over 
the same intermediate wavelengths. During this phase the magnetic energy exponentiates rapidly, with an e-folding time controlled by 
shearing at the smallest available scales. After saturation, the magnetic energy spectrum conforms to the kinetic energy spectrum, such that 
Elf (&) ~ 4Ex(k) over the intermediate wavelengths. The kinetic energy spectrum after the dynamo completion has a slightly shallower 
slope than during the kinematic phase, but is still consistent with a 5/3 spectral index. 



(Kazantsev 



1968). This model predicts that the power 



spectrum oi magnetic energy peaks at the resistive scale 
£„ and obeys a power law k 3 / 2 at longer wavelengths. The 
kinematic phase ends when the magnetic field acquires 
sufficient tension to modify the hydrodynamic motions, 
after which time a dynamical balance between kinetic 
and magnetic energy is established. 

Numerical simulations of MHD turbulence are typi- 
cally limited to magnetic Prandtl numbers Pm = v jr\ w 
1. However, neutron star material is characterized by 
Pm 3> 1, with the viscous cutoff due to neutrino diffu- 
sion occurring at aro und 10 cm, while the resistive scale is 
significantly smaller ( Thompson fc Duncanfll993| ) . How- 
ever, the disparity between true and simulated magnetic 
Prandlt number does not influence our conclusions. This 
is because dynamos are generically easier to establish in 
the high Pm regime th an the small ( Haugen et al-pOO? 
Ponty fc PlunT^[20lT| ) 



We use an initially uniform, pulsar-level (10 11 G) seed 
magnetic field. This field is sub-dominant to the kinetic 
energy by 10 orders of magnitude, so that the initial 
field amplification is expected to be well-described by 
Kazentsev's model. Indeed, we find that during this 
phase the power spectrum of magnetic energy follows 
fc 3 / 2 (Fig. 2|, peaking at around 10 grid zones, which 




large-eddy turnovers 



Figure 3. Time history of the magnetic energy for a represen- 
tative run at 128 3 , together with the empirical model (Equation 

twith best-fit parameters. The horizontal dashed line indicates 
e magnetic energy, E S at at the dynamo completion. From left 
to right, the vertical dashed lines mark the end of the startup, 
kinematic, and saturation phases. 
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we identify as the effective scale of resistivity. The satu- 
ration process begins at ever-earlier times with increas- 
ing numerical resolution. This reflects the fact that dur- 
ing the kinematic phase, magnetic energy exponentiates 
on a time scale controlled by shearing at the smallest 
scales. In numerically converged runs, full saturation 
occurs with Em ~ O.QEk and is characterized by scale- 
by-scale super-equipartition, with Ek ~ ^Em at all but 
the largest scale. 

3.2. Numerical convergence 

The same driven turbulence model was run through 
magnetic saturation at resolutions 16 3 , 32 3 , 64 3 , 128 3 , 
256 3 , and 512 3 . Another model at 1024 3 was run through 
the end of the kinematic phase, but further evolution was 
computationally prohibitive. Fig. [I] shows the develop- 
ment of Bums, Em, and Ek as a function of time at 
each resolution. 

We find that sufficiently resolved runs (> 512 3 ) attain 
mean magnetic field strengths of 10 16 G within two large 
eddy rotations. All models with resolutions > 32 3 even- 
tually attain mean fields of > 3 x 10 16 G. The saturated 
field strength increases until resolution 256 3 . We find 
that the kinematic growth rate is higher at each higher 
resolution, while the time-scale for the non-linear satura- 
tion converges at 256 3 to roughly five large-eddy rotation 
times, or about 0.5ms for the physical parameters of bi- 
nary neutron star mergers. 

In order to quantitatively describe the time develop- 
ment of magnetic energy EM{t), we describe it with an 
empirical model, 



E M (t) = 



E M (0)e t2 ^ 
Em^-^ 
E n \ + {E sat — E n \)(l 



< t < t x 

h<t< t nl 

_ (3) 

where the 6 parameters (summarized in Tab le [l ) are 
obtained by a least-squares optimization. Fig. [3[shows 
the empirical model given in Equation [3] applied to a 
representative run at 128 3 . The first phase, t < ti is a 
startup transient, and lasts until the hydrodynamic cas- 
cade is fully developed at t\ ~ ieddy The kinematic 
dynamo phase lasts between t\ and t a i, during which the 
magnetic energy exponentiates on the time scale n. At 
t n i, the smallest scales reach kinetic equipartition and the 
growth rate slows. In the final stage, Em asymptotically 
approaches £J sa t on the time-scale t%. We define the dy- 
namo completion time t sat as t n \ + r%. Fig. [4] shows the 
best-fit £? S at) tsat) and t± as a function of the resolution. 

The magnetic energy E SB ± at dynamo completion shows 
signs of converging to a value of 0.6Ek by resolution 
512 3 . The time scale on which the magnetic energy 
asymptotically approaches E sa- t is consistently 3i c ddy 
at different resolutions. The dynamo completion time 
t sa t is numerically converged at w 5t c ddy by 256 3 . The 
best-fit kinematic growth time follows a power law in the 
resolution, n cx N~ - 6 . This is consistent with the value 
of —2/3 expected if the dynamo time is controlled by 
shearing at the smallest scale, the cascade is Kolmogorov 
(i.e. Uf oc £ 1 ^ 3 ), and the viscous cutoff l v occurs at a 
fixed number of grid zones. In that case, t± ~ t u oc 




64 J 128' 1 256'' 

resolution 



Figure 4. Top: Convergence study of the kinematic dynamo 
growth time t\ (blue) and the dynamo completion time t sa t (green) 
defined as t n \ + tj. Bottom: Convergence study of the best- 
fit model parameter £ sa t expressed as the ratio of magnetic to 
kinetic energy Ef^/En- The converged value of the volume- 
averaged Em/Ek ~ 0.6. Nevertheless, at intermediate wave- 
lengths Ej^(K)/Eji(k) ft: 4. As shown in Figure [5] the largest 
scales remain kinetically dominated. This indicates the suppres- 
sion of coherent magnetic structure formation near the integral 
scale of turbulence. 

3.3. Power spectrum of kinetic and magnetic energy 

The time development of kinetic and magnetic energy 
power spectra has been studied for a single run with 
resolution 512 3 . We present three-dimensional, spher- 
ically integrated power spectra with the dimensions of 
ergs/cm 3 /m _1 , defined as 



P K {h) = 



1 

Afc~ 



\f*WVp/2 



keAfc; 



P M {h 



Ak. 



keAfci 



(4a) 
(4b) 



,2/3 



TV- 2 / 3 . 



where the Newtonian versions of kinetic and magnetic 
energy are appropriate since the conditions are only 
mildly relativistic. The definitions in Equations [4] sat- 
isfy jP(k)dk = (E) for P K and P M . Figure |] shows 
the power spectrum of kinetic and magnetic energy at 
various times throughout the growth and saturation of 
magnetic field. During the kinematic phase, the kinetic 
energy has a power spectrum Pxik) °c k~ 5 / 3 consistent 
with the Kolmogorov theory for incompressible hydro- 
dynamical turbulence, while Pjy-(fe, t) oc e*/ ri fc 3 / 2 consis- 
tent with Kazenstev's model. Pm(&) maintains the same 
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Table 1 

Empirical model for magnetic energy growth 



Fit parameter 


Description 


Numerically converged value 


TO 


Startup time-scale 


none, artifact 


fl 


Hydrodynamic cascade fully developed 


^cddy 


n 


Kinematic growth time-scale 


none, oc iV -2 / 3 




End of kinematic phase 


2^cddy 


T2 


Non-linear saturation time-scale 


*nl + T"2 ~ 5t eddy , 256 3 


E S a.t 


Saturated magnetic energy 


0.6E K , 512 3 



Note. — Model parameters characterizing the growth of magnetic energy before and during the dynamo saturation. Numerical 
convergence is attained for the fully saturated magnetic energy E Ba ,t and the dynamo completion time t sa t defined as t n i + T2. 



shape, but exponentiates in amplitude at the time scale 
Ti which is controlled by shearing at the resistive scale. 
According to Kazentsev's model, Pnik) should peak at 
the resistive scale. This is consistent with the observed 
peak in the magnetic energy at roughly 10 grid zones, 
the same scale at which we observe the viscous cutoff. 
This is also consistent with Pm = 1 expected from the 
numerical scheme employed. 

When the magnetic energy at the resistive scale sur- 
passes the level of the kinetic energy at that scale, Pm (fc) 
changes shape. The equipartition scale £k,m moves into 
the inertial range, and migrates to larger scale until full 
saturation occurs with Ik.m ~ L/A. The movement 
of (k,m to larger scale is associated with the formation 
of coherent and dynamically substantial magnetic struc- 
tures of increasing size. The time-dependence has b een 
suggested to be £k,m oc i 3 / 2 ( Beresnyak et al.||2009 ). In 
the fully saturated state, the magnetic held is in scale-by- 
scale super-kinetic equipartition throughout the inertial 
range, with Pm(&) ~ 4Pff(fc). The largest scales remain 
kinetically dominated so that the numerically converged 
saturation level is Em ~ O.QEk- 

4. DISCUSSION 

In this Letter we have determined the time scale and 
saturation level of the small-scale turbulent dynamo op- 
erating in the conditions of binary neutron star merg- 
ers. We have presented numerically converged simula- 
tions showing that magnetic fields are amplified to the 
~ 10 16 G level within a small fraction of the merger dy- 
namical time (t sa t ^5 Emerge) j independently of the seed 
field strength. If hydrodynamical instabilities create fluc- 
tuating velocities on the order of 0.1c as indicated by 
global simulations, then each 1km 3 turbulent volume dis- 
sipates > 10 46 erg of magnetic energy per 0.1 ms. If 
/turb represents the fraction of the merger remnant that 
contains such turbulence, the magnetic energy dissipated 
during the merger is at least 



E, 



diss 



10 49 



erg x 



turb 



io- 



lms 



(5) 



A fraction of that dissipation will occur through mag- 
netic reconnection in optically thin surface layers, sup- 
plying relativistic electrons which synchrotron radiate 
in the merger remnant magnetosphere. If 5% of that 
magnetic energy dissipation creates radiation in the 15- 
150 keV band, then the fluence at 200 Mpc would be 
10 _7 erg/cm 2 , potentially rendering most merging neu- 
tron stars in the advanced LIGO and Virgo detection vol- 
umes detectable by Swift BAT. If so, then merging neu- 



tron stars are accompanied by a prompt electromagnetic 
counterpart, independently of whether a later merger 
phase yields a collimated outflow capable of powering 
a short gamma-ray burst. 

We suggest that merger flares may be present in the 
current sample of short GRBs and may be roughly 
isotropic on the sky since they are seen to distances where 
the cosmological matter distribution becomes homoge- 
neous. Searches for merger flares should seek to identify 
short flares, not unlike soft-gamma repeaters, among the 
short burst population. If mergers also produce short 
GRBs short-hard GRBs, then merger flares may consti- 
tute a precursor component of the emission. 

The presence of strong magnetic fields may also aid 
in the ejection of neutron-rich material from surface 
layers of the merger remnant, possibly enhancing the 
enrichment of inter-stellar medium by r-process nuclei 
( Freiburghaus et al. 1999| Rosswog et al. 19991. En- 
hanced production of r-process nuclei also increases the 
likelihood of EM detection by radio-active decay powered 
afterglows, or "kilonovae" ( Li & Paczynski 1998 Metzger 
|et al.||2010[ ). 

Finally, it has been shown that magnetic fields will 
significantly influence the gravitational wave signature 
and remnant disk ma ss, if they exist at the 10 17 G level 
(Etienne et al. 2012). Such strong fields are unlikely 
in older neutron star binaries, but our results suggest 
they may be revived, albeit at small scales, during the 
merger. To have significant influence, those fields would 
have to fill a considerable fraction of the merger volume. 
As we have shown in this Letter, the overall magnetic 
energy budget is controlled by the prevalence (/turb) and 
vigor (ek ) of the turbulent volumes. This fact motivates 
the use of higher resolution global simulations aimed at 
measuring /turb and £k- 
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